Chemical abundance ant icor relations in globular cluster stars: 
The effect on cluster integrated spectra 
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ABSTRACT 

It is widely accepted that individual Galactic globular clusters harbor two coeval generations 
of stars, the first one born with the 'standard' a-enhanced metal mixture observed in field Halo 
objects, the second one characterized by an anticorrelated CN-ONa abundance pattern overim- 
posed on the first generation, a-enhanced metal mixture. We have investigated with appropriate 
stellar population synthesis models how this second generation of stars affects the integrated 
spectrum of a typical metal rich Galactic globular cluster, like 47Tuc, focusing our analysis on 
the widely used Lick-type indices. We find that the only indices appreciably affected by the 
abundance anticorrelations are Ca4227, G4300, CNi, CN2 and NaD. The age-sensitive Balmer 
line, Fe line and the [MgFe] indices widely used to determine age, Fe and total metallicity of 
extragalactic systems are largely insensitive to the second generation population. Enhanced He 
in second generation stars affects also the Balmer line indices of the integrated spectra, through 
the change of the turn off temperature and - in the assumption that the mass loss history of both 
stellar generations is the same - the horizontal branch morphology of the underlying isochrones. 

Subject headings: galaxies: star clusters: general — globular clusters: general — stars: abundances 



1. Introduction 

A large body of spectroscopic data published 
in the last 10 years has conclusively established 
the existence of primordial surface chemical abun- 
dance variations of C, N, O, Na - and sometimes of 
Mg and Al - in individual Galactic globular clus- 
ters (GCs) of all metallicities, from very metal rich 
objects like NGC6441, to the most metal poor ones 
like M15 (see, e.g., Carretta et al. 2010). When 
considering luminosities below approximately the 
red giant branch (RGB) bump level, within a given 
GC the chemical abundance patterns appear to 
be roughly the same all along the RGB down to 
the turn-off (TO) region (see, i.e. Gratton et al. 
2001; Cohen & Melendez 2005), and show CN 
and ONa (and sometimes also MgAl) anticorre- 



lations for about 2/3 of the stars, superimposed 
onto a standard (i.e. consistent with the distribu- 
tion found in field halo stars) a-enhanced heavy- 
element distribution ([a/Fe]^ 0.3-0.4). The anti- 
correlations display a range of negative variations 
of C and O accompanied by increased N and Na 
abundances, whose extension varies from cluster 
to cluster (Carretta et al. 2005, 2009). The re- 
maining 1/3 of cluster stars shows an homoge- 
neous standard a-enhanced metal mixture. At 
luminosities brighter than the RGB bump, obser- 
vations point to the efficiency of some deep mix- 
ing process that decreases the surface C and in- 
creases N in all cluster stars (Yong et al. 2003; 
Gratton et al. 2004; Sneden et al. 2004), similarly 
to the case of field Halo stars (Gratton et al. 2000), 
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whereas the ONa pattern is essentially unchanged 
compared to lower luminosities. This effect varies 
from cluster to cluster, with a cluster like Ml 3 de- 
pleting C much more efficiently than others (Smith 
& Martell 2003; Sneden et al. 2004). 

The accepted working scenario to explain these 
CNONa anticorrelations prescribes that the stars 
currently evolving in a GC were born with the 
observed CNONa patterns. Intermediate-mass 
asymptotic giant branch (AGB) stars in the range 
~ 3 — 8M Q (and/or the slightly more massive 
super-AGB stars) and massive rotating stars are 
considered viable sources of the necessary heavy- 
element pollution (see, e.g., Ventura & D'Antona 
2005; Decressin et al. 2007). Slow winds from 
the envelopes of AGB (or super-AGB stars) or 
from the equatorial disk formed around fast ro- 
tating massive stars inject matter into the intra- 
GC medium after a time of order 10 6 — 10 s yr, 
depending on the polluter. Provided that a signif- 
icant fraction of the material is not lost from the 
cluster, new stars (second generation, but essen- 
tially coeval with the polluters' progenitors, given 
their short evolutionary timescales) may be able to 
form directly out of pristine gas polluted to vary- 
ing degrees by these ejecta, that will show the ob- 
served anticorrelation patterns. Another byprod- 
uct of this pollution may possibly be an enhanced 
initial He-abundance for these second generation 
stars. Recent numerical models by D'Ercole et al. 
(2008, 2010) have started to explore in quantita- 
tive details this broad picture. 

If Galactic GCs are the local counterpart of 
extragalactic globulars, one should expect similar 
abundance anticorrelations to be imprinted in the 
integrated spectra of GCs populating the haloes of 
external galaxies. Given that GCs are employed as 
tracers of the star formation histories of spheroids 
(see, e.g. the review by Brodie & Strader 2006), 
and their spectra are usually modelled consider- 
ing a simple stellar population (SSP) with a single 
initial chemical composition (and age), it is im- 
portant to assess to what extent these abundance 
anticorrelations modify our intepretation of GC 
integrated spectra when compared to the case of 
a SSP. Also, integrated spectra of Galactic GCs 
are often used to test SSP models (see, i.e., Lee & 
Worthcy 2005; Pcrcival et al. 2009; Walcher et al. 
2009; Vazdekis et al. 2010; Thomas et al. 2011) and 
the interpretation of these comparisons is clearly 



dependent on the effect of these abundance anti- 
correlations on the cluster spectra. 

We address for the first time these questions, fo- 
cusing mainly on the predicted Lick-type absorp- 
tion feature indices, that are routinely employed to 
determine ages and chemical compositions of ex- 
tragalactic stellar populations (see, e.g., Burstein 
et al. 1984; Worthey 1994; Trager et al. 1998; 
Proctor et al. 2004; Puzia et al. 2005; Thomas et 
al. 2005; Graves & Schiavon 2008, and references 
therein). The next section presents our methodol- 
ogy and theoretical modelling, and is followed by 
an analysis of the results and a discussion. 

2. Models 

We have started our analysis by considering two 
en-enhanced SSPs - metal mixture with [a/Fe]~ 
0.4 - with ages t=12 and 14 Gyr, [Fe/H]=-0.7 
(Z=0.008) and initial He mass fraction Y"=0.256, 
as representative of the first generation chemical 
composition in typical metal rich Galactic GCs, 
like 47 Tuc, whose integrated spectrum is often 
used for testing stellar population synthesis mod- 
els. For both ages we have then considered a sec- 
ond generation population whose metal compo- 
sition has C decreased by 0.30 dex, N increased 
by 1.20 dex, O decreased by 0.45 dex and Na in- 
creased by 0.60 dex with respect to the first gener- 
ation a-enhanccd mixture, all other metal abun- 
dances being unchanged. This pattern is typical of 
values close to the upper end of the observed an- 
ticorrelation patterns in Galactic GCs (Carretta 
et al. 2005, 2009). The metal distribution of this 
second generation coeval population has the same 
C+N+O sum and the same Fe abundance (as a 
consequence also the total metallicity Z will be 
practically the same) as the first generation com- 
position, in agreement with spectroscopic mea- 
surements on second generation stars within in- 
dividual Galactic GC 1 . 

For both the 12 and 14 Gyr second generation 
populations we have accounted for two alternative 
values of Y, e.g. F=0.256 - as in the first stel- 
lar generation - and F=0.300, to include a possi- 



1 A well known exception to the constancy of Fe within a 
Galactic GC is u> Cen, that displays a large range of Fe 
abundances in addition to He enhancements and CNONa 
anticorrelations (see, e.g., Bellini et al. 2010, and reference 
therein). 
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blc enhancement of He in second generation stars. 
This latter choice is consistent with constraints on 
the typical enhancement of He in Galactic GCs, 
as determined by Bragaglia et al. (2010). We 
stress that these choices for the chemical compo- 
sition of second generation populations are very 
general, and should give us a realistic estimate of 
their impact on GC analyses using SSPs. The ac- 
curate modelling of individual clusters would re- 
quire, however, much more specific prescriptions, 
given that the distribution of CNONa abundances 
in second generation stars varies on a cluster-to- 
cluster basis. One would need to implement the ef- 
fect of a continuous variations of CNONa elements 
up to the maximum values of the anticorrelations 

- that is cluster dependent - plus the exact num- 
ber distribution of stars along the anticorrelation 
pattern, also cluster dependent. Even this would 
not help a detailed comparison with extragalac- 
tic GCs, given that at present we do not know - 
nor we can predict - the distribution of the abun- 
dance anomalies in these objects. The effect of the 
deep mixing along the upper RGB that modifies 
the CN abundances (that to some extent seems, 
again, cluster-dependent) is not accounted for in 
our calculations (nor in any other existing calcula- 
tions of integrated spectra of GCs), but its effect 
goes in the same direction of the CN primordial 
anticorrelation, for it tends to amplify the range 
of C depletions and N enhancements. 

Overall, our selection of representative extreme 
values of the CNONa variations will provide a 
first important indication of the maximum effect 
of these abundance anomalies on GC integrated 
spectra, and will serve as a guideline to interpret 
the abundance pattern derived from fitting SSPs 
to the observed spectra of Galactic and extragalac- 
tic GCs. We have considered in this first investi- 
gation on the subject one single [Fe/H] value, typi- 
cal of the metal rich subpopulations of Milky Way 

- and of 47Tuc, whose integrated spectrum is a 
benchmark for population synthesis models. The 
sensitivity of the spectra to the abundance anoma- 
lies may be quantitatively different when consid- 
ering different [Fe/H] values. 

For each of the first generation SSPs we have 
then calculated synthetic integrated spectra as 
follows. The underlying isochrones are the BaSTI 
a-enhanced isochrones (with Reimers mass loss 
parameter 77=0.2) by Pietrinfcrni ct al. (2006), 



with [Fe/H]=-0.7 (Z =0.008). The 12 and 14 Gyr 
isochrones representing the second generation 
populations are the same as the first generation 
a-enhanced ones with the appropriate Y abun- 
dance (taken from the BaSTI database), for we 
have verified with additional stellar model calcu- 
lations that the effect of the two different metal 
mixtures on the isochrones is negligible, provided 
that the Fe abundance and the C+N+O sum are 
unchanged (see also Salaris et al. 2006). 

We have then computed a grid of 57 stellar 
synthetic spectra for both first and second gen- 
eration metal mixtures, that cover the gravity- 
effective temperature parameter space spanned by 
the corresponding 12 and 14 Gyr isochrones. The 
model atmospheres have been calculated using the 
code ATLAS 12 (Kurucz & Avrett 1981; Kurucz 
2005; Castelli 2005). Spectra were then computed 
with the code SYNTHE (Kurucz & Avrett 1981; 
Sbordone et al. 2004), for the wavelength region 
3500 to 6000 A and convolved to a spectral reso- 
lution R = X/SX = 10000. The atomic line list 
adopted in the computations is based on the com- 
pilations by Coelho et al. (2005) and Castelli & 
Hubrig (2004). Lines for the molecules C2, CH, 
CN, CO, H2, MgH, NH, OH, SiH and SiO from 
Kurucz (1993) were included for all stars, and TiO 
lines from Schwenke (1998) were included for stars 
cooler than 4500 K. 

From the appropriate isochrone and grid of syn- 
thetic stellar spectra we have calculated the inte- 
grated spectrum for these first and second gener- 
ation populations, employing the Kroupa (2001) 
initial mass function, and computed the corre- 
sponding values / of the standard Lick-type in- 
dices as defined in Table 2 of Trager et al. (1998), 
plus the H7 and US indices defined in Worthey & 
Ottaviani (1997). Of the 21 indices in Trager et al. 
(1998) we had to neglect TiOi and Ti02, that are 
beyond the upper limit of the wavelength range 
of our spectra. The values of I\ st generation and 
Iind generation have been determined as in Percival 
et al. (2009), using the LECTOR program by A. 
Vazdekis 2 and measured directly on our high reso- 
lution spectra, i.e., they are not transformed onto 
the Lick system, for the community is rapidly mov- 
ing to measure indices at higher resolution than 
the original Lick system, especially since the pub- 



2 available at http://www.iac.es/galeria/vazdekis/index.html 
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lication of models employing the MILES spectral 
library (Sanchez-Blazquez et al. 2006; Vazdekis et 
al. 2010). As a test, we have degraded the resolu- 
tion of our spectra to 2.3 A, the resolution of the 
MILES spectral library and repeated the analysis 
described below. Although the absolute values of 
the indices have a different zero point, the size of 
the differential effects due to varying metal mix- 
tures (and He abundances) are unchanged. 

To assess the impact of second generation stars 
on ages and chemical composition of unresolved 
GCs, we have chosen to compare our predic- 
tions with the values from the grid of a-enhanced 
SSPs by Percival et al. (2009), that spans a wide 
range of ages and metallicities. These spectra 
have been calculated employing the same BaSTI 
a-enhanced isochrones (with Reimers mass loss 
parameter 77=0.2) used in this paper, the same 
Kroupa (2001) initial mass function, and the syn- 
thetic stellar spectral library by Munari et al. 
(2005). These SSP spectra we adopted as the 
source of a 'reference' grid of values I of the in- 
dices for first generation SSPs with [Fe/H] between 
— 1.84 and +0.05, and age t between 1.25 and 
14 Gyr. We found small offsets in the reference 
index values for the 12 and 14 Gyr, [Fe/H]=-0.7, 
a-enhanced SSPs, when compared to our new AT- 
LAS^ calculations, due to differences in the stel- 
lar spectral libraries. To ensure these small off- 
sets do not play any significant role in our anal- 
ysis, we have decided to use our new results in a 
purely differential way. For each of the selected 
ages and Y of second generation stars, we have 
determined from our calculations the differences 

AI=l2nd generation list generation- The differ- 
ences Al have been then added to the correspond- 
ing index values for the 12 and 14 Gyr 'refer- 
ence' (first generation) a-enhanced grid, to repre- 
sent the corresponding coeval second generation. 
These are the index values employed in the anal- 
ysis that follows. 

3. Results and discussion 

In the following we compare the values of I for 
the various indices in first and second generation 
populations with the same age, to quantify how 
Lick-type indices respond to the presence of a co- 
eval second stellar generation with CNONa anti- 
correlations in individual GCs. It is important to 



stress from the outset that all comparisons that 
follow give a quantitative estimate of the maxi- 
mum effect on the Lick-type indices of GCs, for 
we are comparing a pure first generation with a 
pure second generation population with extreme 
values of the CNONa anticorrelations. In a real 
GC the distribution of CNONa abundances ranges 
(most probably because of the effect of dilution 
with pristine gas) from the values of the standard 
first generation a-enhanced mixture, to extreme 
values like the ones used in our modelling (and 
also, the specific properties of the CNONa distri- 
bution vary from cluster to cluster). One expects 
therefore that in real GCs the variations Al will 
attain values somewhere within the range derived 
from our analysis. 

We have first analyzed the case of the same ini- 
tial He-abundance in first and second generation 
stars. We have considered the 18 metal line in- 
dices that can be measured in our spectra, and as- 
sessed how significant are the corresponding differ- 
ences Al, as follows. At a fixed age of 14 Gyr and 
for each index, we have estimated the metallicity 
[Fc/H] of the first generation population needed 
to match the / value of the second generation pop- 
ulation. If A7=0, first and second generation I 
values are the same, and the 'correct' metallicity 
[Fe/H] =—0.7 of the first generation stars is re- 
quired to enforce agreement. A value of Al differ- 
ent from zero requires instead a [Fe/H] different 
from —0.7 dex for the first generation stars, in or- 
der to match the indices of the second generation. 
Table 1 reports for each metal index the value 
of A J and the corresponding difference A [Fe/H] 
between [Fe/H] and the 'correct' first generation 
[Fe/H]. If we consider as significant A [Fe/H] dif- 
ferences above 0.1 dex, we find that the only metal 
indices appreciably affected by the abundance an- 
ticorrelations are Ca4227, G4300, CNi, CN 2 and 
NaD. The first four indices are sensitive to the 
C, N, O abundances, while NaD is sensitive to 
the Na and C abundances (Trager et al. 1998; 
Korn et al. 2005; Prochaska et al. 2005; Schiavon 
2007; Graves & Schiavon 2008). Not surprisingly, 
CNi, CN 2 and NaD are the indices most affected. 
They display a large increase in the second genera- 
tion population, and to calculate the correspond- 
ing A [Fe/H] we had to perform a linear extrap- 
olation beyond the upper boundary of our avail- 
able [Fe/H] range. The values of the Ca4227 and 
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Table 1: Variations of the values of the Lick- type indices for our 14 Gyr old first and second generation 
populations, both with y=0.256, and the corresponding A[Fc/H] differences (see text for details) . 



Index A J A[Fe/H] 





-0. 


.068 






H7F 





.274 






CNi 





.084 


1 


.960 


CN 2 





.087 


1 


.324 


Ca4227 


-0. 


.651 


-0 


.306 


G4300 


-0. 


.572 


-0 


.360 


Fe4383 


-0. 


.073 


-0 


.022 


Ca4455 


0. 


.000 


-0 


.000 


Fe4531 


-0. 


.014 


-0 


.006 


C 2 4668 


-0. 


.115 


-0 


.063 


H/3 


-0. 


.032 






Fe5015 


0. 


.056 





.016 


Mgi 


0. 


.000 





.000 


Mg 2 


-0. 


.004 


-0 


.016 


Mgb 


-0. 


.158 


-0 


.039 


Fe5270 





.004 





.002 


Fe5335 





.010 





.005 


Fe5406 


-0. 


.005 


-0 


.004 


Fe5709 


-0. 


.060 


-0 


.068 


Fe5782 





.012 





.022 


NaD 


1. 


.346 





.733 



Note. — Units are A but for the CNI, CN2 indices, whose values are expressed in magnitudes. 
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G4300 indices display an opposite behavior, de- 
creasing in the second generation population. The 
other metal indices show much less significant vari- 
ations, when judged on the basis of the associated 
A[Fc/H]. All Fe and Mg indices are essentially un- 
affected, and the same is true for Ca4455, at odds 
with the behavior of Ca4227. The value of C 2 4668 
is also very weakly changed, possibly because the 
decrease of both C and O in the second generation 
population affect the index in opposite directions 
and the resulting net variation is very small (see, 
e.g. Trager et al. 1998). 

We have then investigated the effect of the 
abundance anticorrelations on the Balmer line in- 
dices and cluster age estimates, by considering the 
H/3-Fe5406 and H/3-[MgFe] diagrams 3 displayed in 
Fig. 1. These are two powerful diagrams to break 
the age-metallicity degeneracy and allow the esti- 
mate of, respectively, age and Fe-abundance, and 
age and total metallicity Z of SSPs (see, e.g., 
Worthey 1994; Thomas et al. 2005; Lee & Worthey 
2005; Tantalo & Chiosi 2005; Percival et al. 2009, 
and references therein). Figure 1 confirms the ro- 
bustness of these diagnostic diagrams. The posi- 
tion of the second generation populations on these 
diagrams is barely distinguishable from their co- 
eval first generation counterparts. We have dis- 
played in Fig. 1 also the H£p-Fe5406 and H7F- 
Fc5406 diagrams, that essentially confirm the re- 
sults obtained for the H/3 index. The H7F index is 
more affected by the chemical composition of the 
second generation population, and the value of the 
index increases, a behavior opposite to the case of 
H/3 and H<5p. Ages from H7F are younger by <~2- 
3 Gyr in second generation stars, but given that - 
as mentioned at the beginning of this section - in 
a real GC the variation AJ of a generic index is 
expected to be somewhere within the range deter- 
mined from our model spectra, the 'real' bias on 
the H7F ages is probably smaller, at the level of 
wl Gyr. 

To summarize, the presence of a second gen- 
eration of GC stars with unchanged He-content 
affects appreciably only the Ca4227, G4300, CNi, 
CN2 and NaD metal indices. Very importantly, 
the variation of Ca4227 goes in the direction to 
mimic a lower Ca abundance, if this index is used 



3 [MgFc]= v 'Mg6 X < Fe >, with < Fc >=(Fe5270- 
Fe5335)/2. 



as a measure of the Ca content. Our results 
support Lee & Worthey (2005) suggestion that 
the effect of a second generation population with 
CNONa anticorrelations may explain the discrep- 
ancies they find when comparing Ca4227, CNi, 
CN2 and NaD index strength (on the Lick IDS 
system) from their a-enhanced SSP models, with 
data for Galactic and M31 GCs. On the other 
hand, age, Fe-abundance and Z inferred from H/3- 
Fe5406 and H/3-[MgFe] diagrams (or the H(5p and, 
to a slightly smaller extent, the H7F counterpart) 
are confirmed once again to be robust and insensi- 
tive to the chemical abundance pattern of second 
generation stars. 

To close our analysis, we consider the case of 
the 14 Gyr second generation population with en- 
hanced He-abundance (F=0.300), also displayed 
in Fig. 1. The values of all metal indices in this 
second generation population are essentially iden- 
tical to the Y=0.256 case, i.e. they are unaffected 
by the enhancement of He, but the Balmer line 
indices are changed by the increase of He. The 
values of H/3, H<5p and H7F for this population 
are essentially identical to the 12 Gyr second gen- 
eration population with F=0.256. As a numer- 
ical test, we have determined from our new cal- 
culations all I values for the 14 Gyr second gen- 
eration population isochrone with F=0.256, but 
using the y=0.300 spectra, and found that both 
metal and Balmer line indices are identical to the 
values obtained from the spectra with the appro- 
priate y=0.256. This means that the changes in 
H/3, H(*>f and H7F compared to the Y=0.256 case 
are due to differences in the underlying isochrone 
representative of the He-enhanced second genera- 
tion population. The TO and main sequence of the 
underlying isochrone are hotter by ~100 K com- 
pared to the y=0.256 isochrone. Given that the 
TO mass at fixed age is smaller (by about 0.07 M Q 
in our case) in the He-enhanced isochrone, and 
the total mass lost along the RGB is the same 
(0.11 M Q ) because of our assumption of the same 
mass loss law in all populations, the mass evolving 
along the HB is smaller, hence the typical T c ff of 
HB stars is higher, again by ^100 K. These higher 
T c ff values for TO and HB increase the value of 
the Balmer line indices. 

As a conclusion, the presence of a second gen- 
eration of GC stars with enhanced He does not 
modify the results obtained for the metal lines for 
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Fe5406 Fe5406 

Fig. 1. — Several index-index grids for our reference a-enhanced SSPs. Metallicity increases from left to 
right (in the bottom two panels: [Fe/H]=-1.01, -0.70, -0.29, +0.05, corresponding to Z=0.004, 0.008, 
0.01, 0.0198), and age increases from top to bottom (in the bottom two panels: t=3, 6, 8, 10, 14 Gyr). Filled 
squares denote the 12 and 14 Gyr first generation SSPs, open squares the corresponding second generation 
SSPs with F=0.256, whilst crosses denote a 14 Gyr second generation SSP with enhanced He (see text for 
details). 
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the case with Y=0.256, but impacts the Balmer 
line indices, through the T c ff change of TO and 
(assuming a similar amount of mass lost along the 
RGB in both first and second generation stars) 
HB stars. At low- and intermediate GC metal- 
licities a small change in the HB evolving mass - 
induced for example by an enhanced He in second 
generation stars - is able produce very blue HB 
morphologies, that can have a huge impact on the 
cluster ages determined from Balmer line indices 
(see, e.g., De Propris 2000; de Freitas Pacheco & 
Barbuy 2000; Lee et al. 2000; Schiavon et al. 2004; 
Percival & Salaris 2011, and references therein). 
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